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Abstract. We show that the Suslov nonholonomic rigid body problem stud- 
ied in [I0irr3ll26l can be regarded almost everywhere as a generalized Chaply- 
gin system. Furthermore, this provides a new example of a multidimensional 
nonholonomic system which can be reduced to a Hamiltonian form by means 
of Chaplygin reducing multiplier. Since we deal with Chaplygin systems in 
the local sense, the invariant manifolds of the integrable examples are not 
necessary tori. 



1. Introduction 

We start the paper with the definition of nonholonomic Chaplygin systems, their 
reductions and Hamiltonization. 

Chaplygin Systems. Suppose we are given a natural nonholonomic system on the 
n-dimensional Riemannian manifold (iV, x) with local coordinates Xi , Lagrangian 
l(x, x) = h Y] Xijiiij — v(x) and fc-dimensional distribution D C TN describing 
kinematic constraints: a curve x(t) is said to satisfy the constraints if x(t) £ -D x (t) 
for all t. The trajectory of the system x(t) that satisfies the constraints is a solution 
to the Lagrange-d'Alembert equation 

i— 1 x 7 

Assume that N has a principal bundle structure ir : N — ► Q = N / © with respect 
to the left action of a (n — fc)-dimensional Lie group © and D is the collection of 
horizontal spaces of a principal connection. Given a vector X x G T X N, we have the 
decomposition X x = X' x l + X x , where X x £ D Xl X x £ V x . Here V x is tangent space 
to the fiber © • x (the vertical space at x). 

Further, suppose that the Lagrangian / is also ©-invariant, i.e., v is a (5-invariant 
function and © acts by isometrics on Riemannian manifold (iV, x) . Then the con- 
strained Lagrangian l c (x,x) = l(x,x h ) induces a well defined reduced Lagrangian 
L : TQ — > R via identification TQ w D/&. The reduced Lagrangian L is of the 
natural mechanical type as well, the corresponding kinetic energy (metric) and the 
potential energy will be denoted by xn and V respectively. 

As a result, equation is ©-invariant and defines a reduced Lagrange-d'Alembert 
system on the tangent bundle TQ (for the details see |15ll31lo]L After the Legendre 
transformation it can be rewritten as the following first-order dynamical system on 
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T*Q 

3H dH 

oPi dqi 

where q — (qi, . . . , qu) are local coordinates on the base space Q and Pi = dL/dqi, 
i = 1, . . . , k are conjugate momenta. The Hamiltonian H = ^ ^ K^PiPj + V(q r ) is 
the Legendre transformation of the reduced Lagrangian L = | K Dijqiqj — V(q). 
The functions II,; are quadratic in momenta and depend on the curvature of the 
principal connection and the metric x. 

The system (N,l,D,<&) is called a (generalized) Chaplygin system (see Koiller 
|15p. as a generalization of classical Chaplygin systems with Abelian symmetries 

m 

Chaplygin Reducing Multiplier. Let Af(q) be a differentiablc nonvanishing 
function on Q. Then, under the time substitution dr = M{q)dt the following 
commutative diagram holds 



TQ{q,q} TQ{q,q>} 



4 i p - 



■Af 2 xnq' 



T*Q{q,P} T*Q{q,p}. 
Here q' — dq/dr. In the coordinates {q,q'} and {q,p}, L and H take the forms 

L * = \ E J^XDiMti ~ v (i) and = \ E 772"^ + ^ (</)> 

respectively. We look for a factor Af(q) such that after the above time substitution 
the equations © take the form 

(3) q[ = ^—, Pi = -^r-, i = l,...,ft. 

dpi dqi 

That is, they become Hamiltonian with respect to the symplectic form Q — ^2 dpi A 
dqi = Af(q)Q, + PidAf A dqi, where f2 = dPi A dqi is the canonical symplectic 
form on T*Q. In nonholonomic mechanics the factor Af is known as the reducing 
multiplier. 

It appears that non-existence of an invariant measure of the reduced system 
is an obstruction to its reducibility to a Hamiltonian form. Namely, suppose that 
the original system @ is transformed to the Hamiltonian form with a reducing 
multiplier JV. Then the system has the invariant measure Af(q) ^Clr (see [221 1111 
E]). According to the celebrated Chaplygin's reducibility theorem (see |S|[71^3 or 
section III. 12 in US]), for k = 2, the above statement can also be inverted: the 
existence of the invariant measure with the density Af(q) implies that in the new 
time dr = JV(q)dt, the system gets the Hamiltonian form 

Necessary and sufficient conditions for the existence of an invariant measure of 
the reduced system are given by Cantrijn, Cortes, de Leon and Martin de Diego 
P3- Recently, a nontrivial example of a nonholonomic system (multidimensional 
generalization of the Veselova problem on noholonomic rigid body motion |25lllUp 
for which the Chaplygin reducibility theorem is applicable for any dimension is 
given by Fedorov and Jovanovic 

Note that there is an alternative (but equivalent) description of the method of 
reducing multiplier. The system can be written as Q, n h{XH , •) = dH(-), where 
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Q n h = f2 + S is a nondegenerate 2-form. Then the Chaplygin multiplier is a function 
N such that the form Afflnh is closed (see [221 El El)- 

Contents of the Paper. In Section 2 we introduce the notion of quasi-Chaplygin 
systems. In Section 3 we give a brief description of the Suslov problem and show 
that it can be considered as a quasi-Chaplygin system. Furthermore, this provides 
an example which can be reduced to a Hamiltonian form via Chaplygin reducing 
multiplier. In this sense, the complete integrability of the reduced Suslov problem 
can be defined in the natural way (Section 4). The Hamiltonian description explains 
the solvability of multidimensional Kharlamova, Klebsh-Tisserand and Lagrange 
cases obtained in The topology of invariant manifolds of the Kharlamova 

and Klebsh-Tisserand cases is studied in Section 5. Finally, we note that the 
Lagrange case can be treated as a generalized Chaplygin system in two different 
ways. In the first approach the reduced system becomes Hamiltonian after the time 
rescaling, while in the second approach the reduced system is already an integrable 
Hamiltonian system, namely the multidimensional spherical pendulum. 

2. Quasi-Chaplygin Systems 

Several nonholonomic mechanical systems have classically been regarded as Chap- 
lygin systems in certain properly chosen local coordinates (an example is the well 
known Chaplygin skate), although, globally, they are not Chaplygin systems in the 
sense of the above definition (e.g., see 

Definition 1. With the above notation, we say that (N, I, D, ©) is a quasi-Chaplygin 
system, if we allow that the sum D x + V x does not span the tangent space T X N on 
some ©-invariant subvariety S C N . 

As for the Chaplygin systems, the Lagrange-d'Alembert equation (JJJ is ©- 
invariant and reduces to the quotient space D/<3 (e.g., see The later has 

a structure of the R fc -vector bundle 

R k — > D/<8 

i 

Q = N/<S 

which is not, in general, diffeomorphic to TQ. Outside N \ S, we can treat the 
system as a usual Chaplygin system that reduces to (D|tv\s)/© ~ T(Q \ (5/(8)). 
The jumping in the rank of the distribution D x + V x leads to several interesting 
properties of the system (see examples given below). 

Locally, the reduced system can be derived by the use of Poincare-Chetayev (or 
Bolzano-Hamel) equations. Consider an open ©-invariant set U C N with local 
coordinates x = (q,g), in which the ©-action is simply a ■ (q,g) = (q,ag), a G ©. 
The Lagrangian I is ©-invariant. Whence l(x,x) = l(q,q, 5 _1 .g). 

Let X\, . . . , X n be linearly independent ©-invariant vector fields on U . Then the 
commutators [Xi, Xj] can be written in the basis Xi, . . . , X n : [Xi,Xj] = J2k c ijXk, 
where structural coefficients = c ij(q) & re ©-invariant functions. Let uj\, .. . ,uj n 
be the quasi-velocities defined by x = (q,g) — ^jWjAj. From the definition of 
u>i, one get the relations qi = Xw=i AijUij, where coefficients Aij = Aij(q) are also 
©-invariant functions. 

Now, write the Lagrangian as a function of q and u>: l(q,uj) = ^q^q^g^g). 
Further, suppose that the distribution D is spanned by Xi, . . . ,Xk- Then, the 
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constraints are 

u k+1 = 0, . ..,w n = 

and {qi, . . . ,qk,uJi, . . . ,ujk} can be regarded as local coordinates on D/<8. With the 
above notation, the Poincare-Chetayev equations of the system (e.g., see |18l I??]): 

5(^)=EHi(^i+^. i = i,...,k, 

together with the kinematic equations qi = Aij(q)u)j, f° rm a closed system 

in variables {qi, . . . , qk, wj., . . . , cj^}. Note that, contrary to Chaplygin systems, the 
matrix (j4y : (o F ))i<i,j<Jt does not need to be invertible for all q. 

3. Suslov Problem as a Quasi-Chaplygin System 

Suslov Problem. Consider the motion of an n-dimensional rigid body around a 
fixed point O in the n-dimcnsional Euclidean vector space (V, (•, •)). Let <E\ t . . . , <£« 
be the orthonormal frame fixed in the body and ei, . . . , t n be the orthonormal frame 
fixed in the space. The configuration space of the system is the Lie group SO{n): 
the element g £ SO(n) maps the moving coordinate system to the fixed one. We 
use the following usual matrix notation. Let 

Ex = (1,0,. ..,0)*, E n = (0, ...,0,1)'. 

We take {(£i, . . . , (£„} for the base of V. Then E±, . . . ,E n and 

^1 (^11; ■ ■ ■ i , . . . , e n — (c n i, ■ ■ ■ , 6nn) j &ij — (^25 ) 

will be the coordinate expressions of (Si, . . . , (£ n and of ci, . . . , z n respectively. The 
matrix g £ SO(n) maps the vectors in the moving frame to the same vectors 
regarded in the fixed frame. Therefore Ei = g ■ a and g = (ex,...,e n ) , i.e., 
9ij = Gij. Note that we can consider the components of the vectors e±, . . . , e„ as 
redundant coordinates on SO(n). 

For a path g(t) £ SO(n), the angular velocity of the body is defined by uj(t) = 
9 1 ' £ so(n). From the conditions — Ei = g ■ ei + g ■ ei, we find that e\, . . . ,e n 
satisfy the Poisson equations 

(4) <=i = —u ■ i=l,...,n. 

The kinetic energy of the rigid body is a left-invariant function on TSO(n) of 
the form ^(Ico,u>), where X : so(n) —> so(n) in non-degenerate inertia operator 
and (•, •) denotes the Killing metric on so(n). For a "physical" rigid body, Tuj has 
the form Ioj + ujI, where / is a symmetric n x n matrix called mass tensor (see 
Further, suppose that the body is placed in a potential field that is invariant 
with respect to the orthogonal transformations which fix e„. In our notation this 
means that the Lagrangian has the form 

(5) l(g,g) = hlcj,ui) -v(e n ). 

The Suslov problem describes the motion of a rigid body with the left-invariant 
constraints (see Fedorov and Kozlov [TT)] ) 



(G) 



{u, Ei A Ej) = 0, 1 < i < j < n - 1. 
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Equivalently, we can say that the velocity g belong to the left-invariant distri- 
bution D g = g -X> C T g SO(n), where 

D = span{^i A E n , £ n _i A E n }. 

The motion of the system is described by the Lagrange-d' Alembert equation Q , 
which in the left trivialization takes the form of the Euler-Poincare-Suslov (EPS) 
equation (see [TU] 1 

(7) (4: (Zw) - [Zu, w] — -77— A e n ,r] ) = 0, for all r, G 0. 
\ ere Ce„ / 

The EPS equation, together with the Poisson equations and the constraints 
© completely describe the motion of the Suslov problem in the variables {ei, . . . , e n , ui}. 
For n — 3 these are the equations of the classical Suslov problem with the following 
nonholonomic constrain: the projection of the angular velocity to the vector 
equals zero (see PSIHEI H]). 

Geometry of the Constraints. Let H = SO(n — 1) be the subgroup of SO(n) 
with the Lie algebra t) = = so(n — 1): 

"V' T). °=<» »>■ 

Through the paper we shall simply write SO(n — 1) instead of H. We have the 
following simple geometrical lemma. 

Lemma 1. The left action of SO(n — 1) represents the rotations of a body around 
the vector t n fixed in the space, while the right action represents the rotations around 
the vector (5 n fixed in the body. 

Consider the left action of SO{n — 1) on SO(n) and the principal bundle 

SO(n-l) — > SO{n) 

(8) 4 

S" 1 - 1 = SO(n)/SO(n-l) 

According to LemmaQ] the sphere S 1 " -1 (the space of cosets SO(n — 1) ■ g) can be 
identified with the vector e n via bijection e„ < — > SO(n — 1) • (ei, . . . , e n ) . 

The distribution D C TSO(n) is invariant with respect to the left SO(n — 1)- 
action on TSO(n), but it cannot be regarded as a connection of the principal bundle 
However D can be seen as a connection almost everywhere on SO(n). 

Theorem 2. The distribution D can be regarded as a principal connection of 
the bundle outside the submanifold {e nn = 0}. 

Proof. Since we deal with the left SO(n — l)-action, the vertical distribution V is 
right invariant: V g — so(n — 1) • g. We shall prove that D g and V g span the tangent 
space at g outside the submanifold {e nn = 0}. 
We have: 

D g + V g ^ T g SO{n) <S=J> + g- 1 ■ so(n - 1) • g ^ so(n) 

g ■ • g^ 1 + so(n — 1) ^ so(n) 

(9) <^=^ g ■ X ■ g^ 1 e so(n — 1), for some X GO. 
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Suppose that there is a vector X € such that g ■ X ■ g 1 6 so(?i — 1). Let 
<? = g~\ Ei = g ■ ei. \i X = x\Ei A £„ + ■■■ + x„_iE n _i A E n , then, by using 
g -1 ■ Ei A Ej ■ g = gj A &j , we get 

g- 1 ■ X ■ g = xi&i A e„ H h x n -ie n -i Ae n = Y Ae n E so(n - 1), 

where Y = xiei + • • • + x n -\e n -\. On the other hand, Y A e„ belongs to so(n — 1) 
if and only if 

(10) Y„ = xiein H h x n _ie n _i )n = and e„„ = 0. 

But e„ n = e n „, hence we proved that © implies e nn = 0. 

To prove the opposite statement, we just note that the relation H10J1 considered 
as an equation in the variables X{ always has a solution. Let us choose xi such 
that l |l OB holds. Then, if e nn = e nn = 0, we conclude that the element X = 
x\Ex A E„ + ■ ■ ■ + x n -±E n -i A E n G d satisfies g^ 1 ■ X ■ g = g ■ X ■ g^ 1 € so(n - 1). 
The theorem is proved. □ 

Reduced Suslov Equation. The Suslov equations Q, 0, ® are SO(n — 1)- 
invariant and equations (JJJ , JHJ together with e n = —u> ■ e n can be viewed as their 
SO(n — l)-reduction to 

S n-i x K «-i{ en;W } pa D/SO(n-l). 

For the simplicity, denote the vector e n by q = (qi, . . . , q n ). Further, suppose 
that 1 preserves the decomposition so(n) = so(n— 1)+5. Note that (so(n), so(n—l)) 
is a symmetric pair, i.e., C so{n — 1). Therefore, if is an eigenspace of X, 

then, in view of the condition pr B [5,Z)] = 0, the EPS equation Q and the Poisson 
equation for e„ become 

dv dv 

(11) (Iu)in = ^-Qn ~ T^—Qi, 

oqi oq n 

n-l 

(12) qi = -Ui n q n , i=l,...,n-l, q n = }^u)j n qi- 

i=l 

Contrary to the general case, the reduced Suslov problem Ijllt , l|12() preserves the 
standard measure in the variables {qi, . . . , q n ,u>i n , t02n, ■ ■ ■ , k>n-i,n}- By the Euler- 
Jacobi theorem (e.g., see PP), for the integrability (more precisely, solvability by 
quadratures) of the reduced problem Qlip. 112|) . we need 2n — 5 additional integrals 
which are independent of the energy integral h(Iu),ui) + v(q). In particular, for 
n = 3 we need only one additional integral (see |14lllKl lT]'). 

However, as it is shown in we do not need 2n — 5 integrals to solve the 
multidimensional variants of Kharlamova, Klebsh-Tisserand and Lagrange cases 
of the Suslov problem. Below we will give the natural definition of a complete 
integrability of , (|12|l such that the examples studied in ^2] provide the simplest 
completely integrable cases. Note that, in general, it is still not clear how to define 
the notion of a complete integrability (in the sense of the Liouville theorem) of 
nonholonomic systems (see I16II2] 1 ). 

Chaplygin Reduction. The Lagrangian function I and the distribution D are 
both invariant with respect to the left SO(n — l)-action on TSO(n). As follows 
from TheoremEl the multidimensional Suslov problem is a quasi-Chaplygin system. 
Furtermore, after the appropriate time rescaling, the system becomes Hamiltonian. 
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The reduced space (SO(n) \ {e nn — 0})/SO(n — 1) is the union of two half- 
spheres with q n > and q n < 0. On the half-spheres we can use coordinates 
q = (qi, . . . , g n -i) within the ball 

B = {q eR"- 1 \ (q,q) = q! + ■ ■ ■ + ql_, < 1}. 

Using JSJ) and l|12[l we can write down the reduced Lagrangian and Hamiltonian 

(13) L± = \ l-\q,q) ^ 4) ~ V± ^ H ± = \{ l -^l)){ AP i P ) + V ±(l)i 



where V±(q) = v(q,q n ), q n = ±i/I - (q,q), the metric J is given by J tj = Z m j n , 
P = -Lj q and A = J" 1 . 

We do not need to find the curvature of the connection to find the reduced 
system. In view of (|11|) . after straightforward computations, we can get: 

Theorem 3. Suppose that the inertia operator X preserves the decomposition so(n) = 
so(n — 1) + 5. After the Chaplygin reduction outside submanifold {e nn = 0}, the 
multidimensional Suslov problem Q), Dp, takes the following form on T*B 

(14) « = l§r' p = -^§ q ± + ( AP 'i) p -( P ' AP )i- 

Under the time substitution dr = q n dt and an appropriate change of momenta 
p = q n P, the reduced system {12$ becomes a Hamiltonian system describing a mo- 
tion of a particle within the ball B, 

, . . dHl A . dHl 8V „ dV 

15 q' = —±=Ap, p> = -—± = -— Jq " = - 

ap oq oq Oq 

where ff±(g,p) = \{Ap,p) + V±(q). 

LL Systems as a Quasi-Chaplygin Systems. We can consider Suslov-type 
problems on other Lie groups as nonholonomic systems with left-invariant La- 
grangians and left-invariant nonintegrable distributions (see [151 112| and references 
therein) . Suppose we are given a natural nonholonomic system with a left-invariant 
Lagrangian I and a left-invariant distribution D on a Lie group G (so called LL sys- 
tem). Let d be the restriction of D to g = Tj^G. Futher, suppose that there is a 
subgroup H of G such that for its Lie algebra \) and c) hold = f) + 9, f) fl 5 = 0. 
Then LL system is a quasi- Chaplygin system (G,l, D, H) with respect to the left 
action of H. The system reduces to the dim O-dimensional vector bundle over the 
homogeneous space G/H. As for the Suslov problem, the subvariety S C G where 
the vertical spaces V g = f) • g and the distribution D g = g ■ do not span the tangent 
spaces T g G is given by the equation g • Q • g~ x + I) Q. 

4. Integrability 
In view of j!2jl. Ijl5[l . we can identify 
(16) -w in < — ► (Ap)i, i = l,...,n-l. 

and write (|12[> in the form 



<9v <9v 
' dq dq n 

(17) ? = <?„Ap, q n = -(Ap, q). 



oq dq„ 
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The phase space M = 5" _1 {q} x K" 
has the following natural decomposition 

M = S n ~ x {q} x 

q 



L {p} of the reduced Suslov problem fTTjl 



n-1 {p} » T*BU A T*B, 



where A be the boundary oiT*B. 

According to Theorem in the new time r and outside the domain A the 
dynamics is Hamiltonian. Note that for q n > 0, time r has the same direction 
as t, while for q n < the direction is opposite. Therefore, the closure of flows of 
Hamiltonian vector fields Xh* arid Xh* to T*B recovers the dynamic of (|17|l in 
the following manner. Let 



(q,p) G A\(Ap,q) = 0, 



dv 







q=(9,0) 



be the set of equilibria points of the system Ijl7|l with q n = 0. 

Let (<zo,Po) G A\S and let (g + (r),p+(r)) and (g_(r),p_(r)) be the trajectories 
of Hamiltonian vector fields Xjj* and Xh* such that 



lim (g + (r),p+(r)) 

T— > + T 

Then, invrting the quadrature 



lim (g_(r),p_(r)) 

•— T 



(go,po)- 



r V 1 ~ (9±( T ):9±( T )) 



t-tn 



we find f = t(r) and the dynamics of q = (g, g„) in the original time t. 

The points (<?,f>) G S that correspond to equilibrium points of H17f) are reached 
in an infinite time t. They not need to be the equilibrium points of vector fields 
Xh* and Xh * • 

Definition 2. We shall say that the reduced Suslov problem l|17|) . is completely 
integrable if one can find n — 1 independent smooth integrals /j : .M — * R, which 
after substitutions q n = ±yl — (g, g) Poisson commute between themselves. 





Figure 1. (a) Complete integrability 
In particular, we can take /i = ^(Ap,p) ^ 



(b) Kharlamova case 
w(q). In the case of integrability, 



due to the Liouville theorem, one can (locally) find trajectories (9±(t),p±(t)) by 
quadratures (e.g., see P). On the other hand, the topological part of the theorem 
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can not be directly applied. Namely, if a compact connected component C c of the 
regular invariant set 

(18) M c = {(q,p)eM|/i(q,p)=c i) i=l,...,»-l} 

does not intersect A, then C c is an (n — l)-dimcnsional torus with a uniform quasi- 
periodic dynamics in r and therefore with a non-uniform quasi-periodic motion in 
the original time t. However, if C c intersect A, it may have a quite complicate 
topology (as an illustration, see Figure la). In the three-dimensional case, the 
topological structure of invariant manifolds of several integrable variants of the 
Suslov problem was studied by Tatarinov [21] and Okuneva |19l |20j ■ 

Reconstruction. To reconstruct the motion (g(t),g(t)) on the whole phase space 
D, we have to solve the Poisson equations (@J) for ei,...,e n _i, i.e., to find all 
trajectories in D that projects to the given trajectory (e n (t),u(t)) — (q(i), Ap(t)). 
If (e n (to), u) (to)) is an equilibrium point or if (e n (t),id(t)) is a periodic orbit, then 
the invariant set ir~ 1 {(e n (t) J ui(t))} C D is foliated with invariant tori of maximal 
dimension rank SO(n — 1) or rank SO(n — 1) + 1, respectively (e.g., see Zenkov and 
Bloch |23 and references therein). 

5. Examples of Topology of Invariant Manifolds 

Systems with the Time Symmetry. Suppose that the potential satisfy the 
condition V+(q) — V-(q) — V(q). Then the reduced Suslov equations l|ll|) . (|12H . 
i.e., H17[). are invariant with respect to the transformation 1 i— ► —t, q n — > —q n , and 
£ is simply given by 

£ = {(p,g)G A | (Ap,q) = 0}. 

Thus, in the time t, we are going along the trajectory (q + (r),p + (r)) of Xh^ 
until we reach A\ E. Then we continue to go along the same trajectory, but in the 
opposite direction. If we reach again the boundary A\ S, we get a closed trajectory. 

As an example, consider the Suslov problem with a rigid body inertial operator 
Xlo = lid + luI, I = diag(ii, . . . , /„) and a quadratic potential w(q) = C\q\ + • • • + 
C n -iq n -i + \(B x q\ + ■■■ + B n ql). Then 

H* = HI = - ( —^—p\ + ■■■ + pi ,) + 

(19) +C m + ■■■ + C„_!5 n _i + i ((B x - B n )ql + ■■■ + (B n _i - B^q 2 ^) . 
The functions 

fi(q,p) = p\ + 2d(Ii + I n )qi + (It + I n )(Bi - B n )qf, i = 1, . . . ,n - 1 

Poisson commute and they also commute with the Hamiltonian (|19fl . Moreover, 
they are smooth functions on M. Thus, the system (fT7|) is completely integrable. 

Furthermore, we can consider the flow of and integrals fa on the whole 
symplectic linear space M n_1 {?} x M n ^ 1 {p}. Let T c = {/, = a} C R 2n ~ 2 be an 
invariant submanifold of Xh*. Then the invariant set M c = {fi = Ci} a M is & 
two-fold covering of 

t; = % n {q\ + ■ ■ ■ + ql_ x < 1} 

determined by the multivalued function q n = ±^/l — (q, q). The branching points 
of the covering correspond to zeros of q n , i.e., to dT* w A n M c - 
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The Kharlamova Case. Let w(q) = C\q\ + ■ ■ ■ + C n -iq n -i (multidimensional 
Kharlamova case). This is a potential for a rigid body placed in a homogeneous 
gravitation field with the mass center orthogonal to (£„. The trajectories of the 
system can be found by quadratures and can be expressed in terms of elliptic 
functions of time t (see Jovanovic |13|L 

A generic invariant set T c is diffeomorphic to R™ -1 and T* is a disjoint union of 
2 l copies of {n — l)-dimensional closed balls (the number I, < I < n — 1 depends 
on the choice of constants Ci). Therefore, connected components of A4 C are spheres. 
Trajectories of the vector field Xh* pass through A. Therefore, apart from a finite 
number of equilibrium points and their asymptotic trajectories, all the trajectories 
of the Suslov problem in the original time t will be closed (Figure lb). Thus we get 

Theorem 4. In the Kharlamova case, the phase space M of the reduced problem is 
almost everywhere foliated with (n— 1) -dimensional spheres. The distribution D is 
filled up with conditionally-periodic trajectories of maximal dimension [n/(n— 1)]+1. 

Klebsh Tisserand Case. For generic values of the constants Ci and Bi the topol- 
ogy of invariant manifolds is much more complicated. Consider the quadratic poten- 
tial v(q) = ^(Biqf + - ■ ■+B n q'^ l ) (multidimensional Klebsh-Tisserand case |13p with 
Bi > B n , i = 1, . . . , n— 1. Then T c is the (n— l)-dimensional torus S Cl x • • • x S Cn _ 1 , 
where S Ci = {fi = Ci} are circles in the planes M. 2 {pi,qi}. (If some of a vanish, 
then the dimension of the tori decreases.) Let fi be angular variables on the circles 
S Ci - 

I C ' 

Pi = ^Cicos^i, q% = \ — sin</? 4 , Xi = (Ii + I n ){Bi - B n ), i = l,...,n-l, 

V % 

Then the Hamilton equations on T c take the form 



fon\ ^ n Bi-B n 

(20) d7 = a = Vl7TXr' • = 

and the subset T* C T c is given by equation 

(21) — sm ipx H 1 sin <Pn-i < 1- 

Using the relation 112 1L we can describe the topological structure of the invariant 
manifolds M. c of the Klebsh-Tisserand case. Namely, the condition 

(22) SL + . . . + f2zl < i 

implies that q n ^ on T c . As a result, the following theorem holds (see |13p. 

Theorem 5. In the domain defined by the condition ]2>ty - invariant manifolds 
are disjoin unions of two l-dimensional tori (I < n — 1) with non-uniform 
conditionally-periodic motions. 

For — + •••+ > 1, the boundary of T* is not empty and there are several 
cases depending of the values of Ci . We quote some of them, 
(i) < ci < H\, . . . ,0 < c„_i < H n -i. 

In this case T* is obtained by removing a disjoint union 
of 2 n_1 open balls Bi from the torus T c . The invariant set M. c is diffeo- 
morphic to 

(r;*- 1 \ UiBi) n Ul9Sl (r; 1 - 1 \ u^). 
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For n = 3 this manifold is a sphere with five handles (see Figure 2). 

(ii) Cl > XI, . . . , C„_i > Xn-l. 

7^.* is a disjoint union of 2™ _1 closed balls. Similarly to the Kharlamova 
case, M. c is diffeomorphic to disjoint union of 2™ _1 spheres S 1 ™ -1 . 
(iii) ci > xi,0 < c 2 < x 2 , . . .,0 < c„_i < x n _i. 



The subset T* is diffeomorphic to two copies of T" 



[0,1]. As a 



result, A4 C is a two-fold covering of 2T" x [0, 1] with branching point on 



the boundaries T' 
dimensional tori T 




{0, 1}, i.e., M. c is a disjoint union of two (n — 1)- 




Mc 



a2 



bi 



b3 



b4 



a2 -|- C2 



bi 



b2 



to 



b4 



ai 



ci 



Figure 2. Klebsh-Tisserand Case 



Analysis of other cases goes along similar lines. It is desirable to have a complete 
description of the bifurcation diagram of the mapping F = (fi, . . . , / n _i) : A4 — > 
R™ -1 , as given for n = 3 in |19) . 

For generic values of the frequencies f2j, the trajectories of l|20|) are dense in T c 
and, in the domain +• • •+ e "~ 1 > 1, a generic trajectory intersects the boundary 
dT* w A n A^ c . Therefore, in the original time t almost all the trajectories are 
closed. Thus, concerning the reconstruction problem, it would be interesting to 
describe qualitative behavior of the motions over invariant sets M. c , where ^- + 
• • • + °™~\ < 1 • A similar problem for the Fedorov-Kozlov integrable case of the 
Suslov problem (the original constraints © are relaxed and v(e n ) = 0) is studied 
in |25| . 

Lagrange Case. Consider the EPS equations Q, (0J, JJjJ describing the mo- 
tion of a dynamically symmetric heavy rigid body. Namely, assume that / = 
diag(Ji, . . . ,Ii,I n ) and that the mass center lies on the symmetry axis (5 n . This 
gives us a nonholonomic version of the multidimensional Lagrange top considered 
by Beljaev The Lagrangian of the system has the form 

(23) l(g,g) = Ii(pr ao ( n _ 1) w.pr^.x) w) + + J re )(pr, w,pr B w) - ee nn . 

The reduced system l|17ll is integrable in view of the Euler-Jacobi theorem and 
in the sense of Definition as well. Namely, the functions fij = qipj — qjPi : 
A4 —> R are integrals of (|17|l . There are In — 5 independent integrals of this 
form, which, together with the energy + I n )(p,p) +eq n , implies that the phase 
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space M. is foliated with two-dimensional invariant manifolds. Since the system 
has an invariant measure, it is integrable by the Euler-Jacobi theorem p^. On the 
other hand, the functions /y, H± ensure the non-commutative integrability (and, 
therefore, the usual complete integrability 0) of the Hamiltonian flows Xhi within 
T*B. 

The geometrical meaning of the integrals /y is that the trajectories q(r) take 
place over invariant two-dimensional planes, depending on the initial conditions. 
After fixing the plane, the problem becomes the usual two-dimensional problem of 
the motion of a particle in a central potential force field. 

Reduction to the Spherical Pendulum. The nonconstrained Lagrange top sys- 
tem is completely integrable (see [3], the Lax pair for the system is given in It 
appears that the Suslov problem Q), l@J), © can be regarded as a subsystem of the 
Lagrange top. Namely, the Lagrangian (|23l) is invariant with respect to both left 
and right actions of the Lie subgroup SO{n — 1). Let us consider the right action. 
(Note that D is a principal connection of the bundle (JHJ) given by the right action 
of SO(n — 1).) The momentum map of the system is <p{g, g) — 2I\ w so (n-i) tJ - This 
gives us the conservation law P r so (ri-i) = 0, In particular, the distribution 
D = cf)~ 1 (0) is an invariant submanifold of the Lagrange top system. 

Therefore, the right SO(n — 1)-Chaplygin reduction of the Suslov problem to 
TS 71 ^ 1 coincides with the Lagrange-Routh reduction of the Lagrange top for zero 
value of the momentum mapping (see pQ, page 87, Theorem 13). According to 
Lemma^ the sphere S*™ -1 = SO(n)/ SO(n— 1) can be identified with the positions 
of the vector (£„ considered in the fixed frame. The reduced Lagrangian has the 
form 

L = -(h + I n )(el„ + ■■■ + e 2 nn ) - ee nn . 

Hence, the reduced system is integrable: it is the multidimensional spherical pen- 
dulum. 
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